MADPH 01-1226 
BNL-HET-01/15 
MSUHEP-10510 
DFTT 13/2001 
hep-ph/0105129| 



H + 2 JETS VIA GLUON FUSION 



V. Del Duca a , W. Kilgore 6 , C. 01eari c , C. Schmidt d and D. Zeppenfeld c 

I.N.F.N., Sezione di Torino via P. Giuria, 1 - 10125 Torino, Italy 
^ Physics Department, Brookhaven National Laboratory, Upton, New York 11973, U.S.A. 

( c ) Department of Physics, University of Wisconsin, Madison, WI 53706, U.S.A. 
(d) Department of Physics and Astronomy, Michigan State University, East Lansing, MI 

48824, U.S.A. 



Abstract 



Real emission corrections to gg — > H, which lead to H + 2 jet events, are calculated 
at order a^. Contributions include top-quark triangles, boxes and pentagon diagrams 
and are evaluated analytically for arbitrary top mass m(. This new source of H + 2 jet 
events is compared to the weak-boson fusion cross section for a range of Higgs boson 
masses. The heavy top-mass approximation appears to work well for intermediate 
Higgs-boson masses, provided that the transverse momenta of the final-state partons 
are smaller than the top-quark mass. 



1 Introduction 



Gluon fusion via a top-quark loop is expected to be the most copious source of Higgs 
bosons in high energy hadronic collisions. Because of its importance for Higgs searches at 
the CERN LHC [p], §], precise knowledge of higher order QCD corrections will be crucial to 
predict event characteristics for a Higgs signal, and to extract Higgs properties, like Higgs 
couplings to gauge bosons and fermions jj. 

For the measurement of Higgs boson couplings, the separation of different Higgs produc- 
tion processes is necessary. Weak-boson fusion (WBF), i.e. the radiation of a Higgs boson 
from a t-channel W or Z boson in electroweak quark-quark scattering, is characterized by the 
presence of two forward jets. Such H + 2 jet events can also arise from 0(a 2 s ) corrections to 
gluon fusion. For a measurement of WBF cross sections, the "background" from this latter 
process must be known, i.e. we need to calculate the cross section and event characteristics 
of H + 2 jet events arising from gluon fusion. 

Next-to- leading order (NLO) QCD corrections to the inclusive gluon-fusion cross section 
are known to be large, leading to a X-factor close to two [|J. Because the lowest order 
process is loop induced, a full NNLO calculation would entail a three-loop evaluation, which 
presently is not feasible. Fortunately, for the intermediate Higgs mass range, which is the 
focus of present interest, the Higgs boson mass vriu is small compared to the top-quark pair 
threshold and the large m t limit promises to be an adequate approximation. Consequently, 
present efforts on a NNLO calculation of the inclusive gluon-fusion cross section concentrate 
on the mj— >oo limit, in which the task reduces to an effective two-loop calculation ||. In 
order to assess the validity of this approximation, gluon-fusion cross-section calculations, 
which include all finite m t corrections, are needed. Of particular interest are phase space 
regions where one or several of the kinematical invariants are of the order of, or exceed, 
the top-quark mass, i.e. regions of large Higgs boson or jet transverse momenta, or regions 
where dijet invariant masses become large. 

In this letter we present first results of such a calculation, namely the evaluation of the 
real emission corrections to gluon fusion which lead to if + 2 parton final states, at order 
«g. This includes the subpro cesses 

qQ -> qQH , qg -> qgH , gg -► ggH , (1.1) 

and all crossing-related processes. The corresponding scattering amplitudes, involving tri- 
angle, box and pentagon diagrams with a top-quark loop, are calculated analytically for 
arbitrary values of the top-quark mass. 
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Figure 1: Samples of Feynman graphs contributing to H + 2 jet production via 
gluon fusion. 

2 Calculation 

In Fig. |l] we have collected a few representative Feynman diagrams that occur in Higgs 
production plus two jets, at order . In our calculation, the top quark is treated as massive 
but we neglect all other quark masses, so that the Higgs couples to gluons only via a top- 
quark loop. Typically we have a ggH coupling through a triangle loop (Fig. |l] (a)), a 
gggH coupling through a box loop (Fig. [I] (b)) and a ggggH one through a pentagon loop 
(Fig. 0(c)). 

The number and type of diagrams can be easily traced back if one thinks to insert these 
Higgs-gluon "vertices" into the tree-level diagrams for 2^2 QCD parton scattering. 
In the following counting, we exploit Furry's theorem, i.e. we are counting as one the two 
charge-conjugation related diagrams where the loop momentum is running clockwise and 
counter-clockwise. This halves the number of diagrams. In addition, the crossed processes 
are not listed as extra diagrams, but are included in the final results. 

1. qQ — > qQH There is only one diagram obtained from the insertion of a triangle loop 
into the tree- level diagram qQ^qQ, i.e. Fig. |l] (a). 

2. qg — ► qgH At tree level, there are 3 diagrams contributing to the process qg qg: 
one with a three-gluon vertex and two Compton-like ones. Inserting a triangle loop 
into every gluon line, we have a total of 7 different diagrams. In addition, we can insert 
a box loop into the diagram with the three-gluon vertex, in 3 different ways: the 3! 
permutations of the 3 gluons are reduced to 3 graphs by using Furry's theorem. In 
total we have 10 different diagrams for the qg — > qgH contribution. 

3- gg ggH Four diagrams contribute to the tree-level scattering process gg—>gg: a 
four-gluon vertex diagram and 3 diagrams with two three-gluon vertices each. Inserting 
a triangle loop in any of the gluonic legs, gives rise to 19 different diagrams. The 
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insertion of the box loop in the 3 diagrams with three-gluon vertices yields another 18 
diagrams. Finally, there are 12 pentagon diagrams (corresponding to 4! permutations 
of the external gluons, divided by 2, according to Furry's theorem). 



The procedure for computing these diagrams is outlined below. A detailed discussion is 
postponed to a later paper ||. 



Tensor integrals Triangle and box integrals are evaluated in terms of Passarino-Veltman 
Cij and Dij functions 0. Gluon polarization vectors are expressed in terms of external 
momenta, which substantially simplifies expressions. Helicity techniques are employed for 
the numerical evaluation of all amplitudes involving external quarks. 

In case of the pentagon diagrams appearing in the ggggH amplitudes, the four external 
gluon-momenta form a basis in terms of which we expand the gluon polarization vectors. 
In this way, all the scalar products appearing in the numerator of the tensor integrals are 
products of the external momenta and the loop momentum (k) . These dot products are then 
written as a difference of two propagators, plus a /^-independent term. In this way, tensor 
pentagon integrals are transformed into a linear combination of tensor boxes and of scalar 
pentagons. 

Scalar integrals All the scalar integrals needed for the calculation are finite in D = 4 
dimensions, due to the presence of the top-quark mass. No further regulator is required. 
Scalar triangles and boxes have been known for a long time in the literature || and efficient 
computational procedures are available |J. Following the procedure outlined in Refs. [ ID . 
we have expressed all scalar pentagons as linear combinations of scalar boxes. 

Checks We were able to perform two different kinds of checks on the analytic amplitudes 
we computed: 

1. Gauge invariance: an important test of the calculation is the gauge invariance of the 
amplitudes for qg —>■ qgH and gg —>■ ggH. This test was performed both analytically 
and numerically, in the final Fortran program. 

2. Large rat limit: the amplitudes for Higgs plus two partons agree in the large mt 
limit with the corresponding amplitudes obtained from the heavy-top effective La- 



grangian [|11]]. This was done numerically by setting ra t = 3 TeV. We found agreement 
within a few percent, when the Higgs mass is varied in the range 100 GeV < ran < 
700 GeV. 
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3 Discussion of the results 



The gluon-fusion processes at O(aj), together with weak-boson fusion (qQ — > qQH pro- 
duction via t-channel exchange of a W or Z), are expected to be the dominant sources of 
H + 2 jet events at the LHC. Cross sections for the former diverge as the final-state partons 
become collinear with one another or with the incident beam directions, or as final-state 
gluons become soft. A minimal set of cuts on the final-state partons, which anticipates 
LHC detector capabilities and jet finding algorithms, is required to define an H + 2 jet cross 
section. Our minimal set of cuts is 



p Tj >20GeV, |ty|<5, 



Rjj > 0.6, 



(3.1; 



where pxj is the transverse momentum of a final state jet and Rjj describes the separation 
of the two partons in the pseudo-rapidity r\ versus azimuthal angle plane. 



In our simulation, we used the CTEQ4L set for parton-distribution functions with 
a factorization scale equal to y/prx pr2 and we fixed a s = 0.12. We postpone any discussion 
on factorization- and renormalization-scale dependence to a further paper |J. 
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Figure 2: H + 2 jet cross sections in pp collisions at ^/s = 14 TeV as a function of 
the Higgs boson mass. Results are shown for gluon-fusion processes induced by a top- 
quark loop with mt = 175 GeV and in the mt — > oo limit, computed using the heavy- 
top effective Lagrangian, and for weak-boson fusion. The two panels correspond to 
two sets of jet cuts: (a) inclusive selection (see Eq. ( \3. 1\ ) ) and (b) WBF selection 
(Eqs. pJl) and pJjj. 



Expected H + 2 jet cross sections at the LHC are shown in Fig. ^ as a function of the 
Higgs boson mass, mn- The three curves compare results for the expected SM gluon-fusion 
cross section at m t = 175 GeV (solid line) with the large mt limit (dotted line), and with 
the WBF cross section (dashed line). Error bars indicate the statistical errors of the Monte 
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Carlo integration. Cross sections correspond to the sum over all Higgs decay modes: finite 
Higgs width effects are included. 



Figure [2] (a) shows cross sections within the minimal cuts of Eq. ( |3.1| ). The gluon-fusion 
contribution dominates because the cuts retain events with jets in the central region, with 
relatively small dijet invariant mass. 

In order to assess background levels for WBF events it is more appropriate to consider 



typical tagging jet selections employed for WBF studies Hl3|l . This is done in Fig. || (b) 
where, in addition to the cuts of Eq. ( |3.1|) , we require 



\Vji ~ Vji\ > 4 - 2 > Vji ■ Vj2 < 0, rrijj > 600 GeV, (3.2) 

i.e. the two tagging jets must be well separated, with three units of pseudo-rapidity be- 
tween the jet definition cones, they must reside in opposite detector hemispheres and they 
must possess a large dijet invariant mass. With these selection cuts the weak-boson fusion 
processes dominate over gluon fusion by about 3/1 for Higgs boson masses in the 100 to 
200 GeV range. This means that a relatively clean separation of weak-boson fusion and 
gluon-fusion processes will be possible at the LHC, in particular when extra central-jet- veto 
techniques are employed to further suppress semi-soft gluon radiation in QCD backgrounds. 
A suppression by a factor three of gluon fusion as compared to WBF cross sections is to be 



expected with a central-jet veto []13" |. 



A conspicuous feature of the H + 2 jet gluon-fusion cross sections in Fig. is the threshold 
enhancement at ran ~ 2m t , an effect which is familiar from the inclusive gluon-fusion cross 
section. Near this "threshold peak" the gluon-fusion cross section rises to equal the WBF 
cross section, even with the selection cuts of Eq. ( |3.2| ). Well below this region, the large 
rrit limit provides an excellent approximation to the total H + 2 jet rate from gluon fusion, 
at least when considering the total Higgs production rate only. Near top-pair threshold the 
large m t limit underestimates the rate by about a factor of 2. 

A somewhat surprising feature of Fig. ^ (b) is the excellent approximation provided by 
the large m t limit at Higgs boson masses below about 200 GeV. Naively one might expect 
the large dijet invariant mass, rrijj > 600 GeV, and the concomitant large parton center- 
of-mass energy to spoil the mt—* oo approximation. This is not the case, however. The 
validity range of the m t — > oo approximation is best appreciated in Fig. ^, where we show 
the transverse-momentum distribution of the harder of the two jets for an intermediate mass 
Higgs boson. At transverse momenta below 200 GeV the large m t limit works extremely 
well. It is the large prj region where this approximation breaks down. Small jet transverse 
momenta but large dijet invariant masses are well modeled by the large m t limit. 
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Figure 3: Transverse-momentum distribution of the hardest jet in H + 2 jet events 



from gluon-fusion processes. Jets are defined via the inclusive selection of Eq. (3.1). 
The two curves are for nit = 175 GeV and for the nit — > oo limit, computed using 
the heavy-top effective Lagrangian. The mass of the Higgs is set to 120 GeV. 



4 Conclusion 

In this letter we have provided first results of a full O(a^) calculation of H + 2 jet cross 
sections, including the analytic top-mass dependence. Event rates are sizable, of order 5 to 
10 pb with minimal jet-selection cuts. With stringent forward-jet tagging cuts, as suggested 
for the Higgs search in weak-boson fusion, rates drop to about a third of the WBF rate, 
which implies that a relatively clean separation of gluon fusion and weak-boson fusion will 
be possible at the LHC. Sizable deviations from the large m t limit are seen for Higgs boson 
masses of order 2mt or for large jet transverse momenta. However, the large m t limit works 
well for parton center of mass energies > 2m t , provided that m# < m t and that the 
transverse momenta remain small compared to m t . 
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